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ABSTRACT 

We present the discovery of fast infrared/X-ray correlated variability in the black- 
hole transient GX 339-4. The source was observed with sub-second time resolution 
simultaneously with VLT/ISAAC and RXTE/PCA in August 2008, during its persis- 
tent low-flux highly variable hard state. The data show a strong correlated variability, 
with the infrared emission lagging the X-ray emission by 100 ms. The short time de- 
lay and the nearly symmetric cross-correlation function, together with the measured 
brightness temperature of ~ 2.5 x 10 6 K, indicate that the bright and highly variable 
infrared emission most likely comes from a jet near the black hole. Under standard 
assumptions about jet physics, the measured time delay can provide us a lower limit 
of T > 2 for the Lorentz factor of the jet. This suggests that jets from stellar-mass 
black holes arc at least mildly relativistic near their launching region. We discuss 
implications for future applications of this technique. 

Key words: black hole physics - stars: winds, outflows - X-rays: binaries - X-rays: 
individual: GX 339-4 



1 INTRODUCTION 

The wealth of multi-wavelength observations of X-ray bina- 
ries (XBs) over the past decade have made clear the ubiq- 
uity of jets in these systems (see Fender 2006 for a review). 
These jets are thought to be the origin of the observed ra- 
dio emission. The soft X-ray flux is generally believed to 
come predominantly from accretion discs around these com- 
pact objects, while the hard X-ray flux is thought to arise 
from a hot Comptonizing corona and/or from the jet itself. 
Recently it has been shown that also the infrared (IR) emis- 
sion includes a substantial contribution from the relativistic 



jet, in the hard states of XBs (e.g. Corbel & Fender 



|Malzac, Merloni fc Fabian 2004] |Russell, Fender fc Jonker| 



2002 



* Based on observations made with ESO Telescopes at the 
Paranal Observatory under programme ID 281.D-5034 
f E-mail: p.casella@soton.ac.uk 



20071. Despite the rapid increase in our phenomenological 



understanding of jets from XBs, we still lack a fundamental 
understanding of how jets are powered and collimated, or 
what the bulk and internal properties of the jets are. 

High-speed simultaneous optical/X-ray photometry of 
three accreting black holes (BHs) opened a new promis- 
ing window. A complex correlated variability in the opti- 
cal and X-ray emission ( |Spruit fc~ Kanbach 2002 j) was seen 
from XTE Jl 118+480, while fast optical photometry of 



SWIFT J1753.5-0127 ( |Durant et al.||2008[ ) and G X 339-4 
(jGandhi et al.|[2008[ ) revealed further complexity. [Malzac 



et al. ( 2004 1 explained the behaviour observed in XTE 



J 1118+480 through coupling of an optically emitting jet and 
an X-ray emitting corona in a common energy reservoir. An 
alternative explanation comes from the magnetically driven 



disc corona model (Merloni, Di Matteo & Fabian 2000) 



magnetic flares happen in an accretion disc corona where 
thermal cyclo-synchrotron emission contributes significantly 
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to the optical emission, while the X-rays are produced by 
Comptonization of the soft photons produced by dissipation 
in the underlying disc and by the synchrotron process itself. 
The two explanations differ substantially in the predictions 
at 1R wavelengths, where a jet appears as the most probable 
origin for the emission (e.g. |Russell et aL|2006| ) . 

The BH candidate GX 339-4 is a recurrent X-ray tran- 
sient ( jMarkert et al.|1973[ ). It has been detected as a highly 
variable source from radio through hard X-rays (see e.g. 
Makishima et al.|[T986] |Corbel et aTlpOOOl |Coriat et aT] 
2009| and references therein). Optical spectroscopy indicates 
a mass function of 5.8 ± 0.5 MQ and a minimum distance 



of 6 kpc ( jHynes et aT1|2003[ 2004). Multiwavelength cam- 
paigns clearly reveal a non-thermal contribution to the in- 
frared (IR) emission in the hard state, most probably arising 



from a compact jet (Corbel & Fender 20021. It is the first 
BH XB for which fast optical/X-ray correlated variability 
was observed i Motch, Ilovaisky fe Chevalier|19"82 1. 

Past variability studies on timescales of several seconds 
have been used to suggest a jet origin for IR emission from 
XBs (e.g. jHynes et al.|2003||Eikenberry et al.|2008| |. In this 
Letter, we report on the first simultaneous fast (i.e. sub- 
second) timing IR/X-ray observations of GX 339-4, during 
its 2008 low luminosity hard state. 



2 OBSERVATIONS 
2.1 Infrared data 

We observed GX 339-4 from ESO's Paranal Observatory 
on 18 August 2008. We obtained fast A's-band photom- 
etry with the Infrared Spectrometer And Array Camera 
(ISAAC) ( jMoorwood et al.| [T998;) mounted on the 8.2-m 
UTl/Antu telescope. The 23"x23" window used encom- 
passed the target, a bright 'reference' star (Ks=9.5) located 
13.6 arcsec south of our target and a fainter 'comparison' 
star (K s =12.8) 8.9 arcsec north-east of GX339-4. 

We used the "FastJitter mode" with a time resolution 
of 62.5 ms. This generated cubes of data with 2500 images 
apiece, and with a 3 s deadtime between cubes. The ULTRA- 
CAM pipelin^Jwas used for the data reduction, after apply- 
ing a barycenter correction for Earth motion. We performed 
fixed-aperture photometry of the three sources (target, ref- 
erence and comparison stars) and used the bright reference 
star for relative photometry of the target and comparison 
stars. The positions of the aperture regions around the tar- 
get and the comparison star were linked to the position of 
the bright reference star to allow for image motion and were 
updated at each time step. The atmospheric conditions were 
good and the resulting light curve for the comparison star 
was consistent with a constant, as expected. By combin- 
ing all 250000 images, we estimate a de-reddened (Ay=3.9, 
A K = 0.114 x ^ v =Q.445- |Cardelli, Clayton fc Mathis|1989 l 
average magnitude of Ks = 12.0±0.2 for GX 339-4, which cor- 
responds to an average flux of F ~ 1.5 x 10 -11 erg s _1 cm~ 2 . 
A sample of the highly variable light curve for GX 339-4 is 
shown in the bottom panel of Fig. [I] 




1 We thank Tom Marsh for the use of the ULTRACAM pipeline 
software i jhttp: / /deneb.astro.warwick.ac.uk/phsaap/software/[ l . 
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Figure 1. Top panel: A sample of the X-ray light curve of GX 
339-4, obtained with the PCA onboard RXTE. The data are 
background subtracted, in the 2-15 keV energy range, at 1-second 
time resolution. Bottom panel: The simultaneous IR light curve, 
obtained with ISAAC. We show the ratio between the source (av- 
erage 4.4 X 10 5 counts/s) and the reference-star (6x 10 6 counts/s) 
count rates in the Kg filter, at 1-second time resolution. The right 
ordinates show the de-reddened flux. We show the typical error 
bars in the top-left corner of each panel. 



2.2 X-ray data 

Simultaneously with the IR observations, GX 339-4 was 
observed with the Proportional Counter Array (PCA) on- 
board the Rossi X-ray Timing Explorer (RXTE). Two pro- 
portional counter units (PCUs) were active during the whole 
observation. The X-ray data span three consecutive satellite 
orbits, for a total exposure of 4.6 ksec. The Binned Mode 
(8 ms time resolution) was used for this analysis, using the 
2-15 keV energy range (channels 0-35). The barycenter cor- 
rection for Earth and satellite motion was applied. Standard 
HEAD AS 6.5.1 tools were used for data reduction. In the 
upper panel of Fig. [l] we show a sample of the light curve, 
corresponding to the second RXTE orbit. Spectral fitting 
with a power-law with photon index 1.6 results in a 2-10 
keV unabsorbed flux of Fx ~ 1.4 x 10~ 10 erg s _1 cm -2 . 



3 CROSS-CORRELATION FUNCTION (CCF) 

Both datasets have an absolute time accuracy better than 
the time resolution used here: ISAAC data have a timing 
accuracy of about 10 ms (the readout time), while RXTE 
data have a timing accuracy of 2.5 /us ( Jahoda et al.|2006 i. 

From Fig. [I] a strong correlation between X-ray and 
IR flux is evident. Both long, smooth variability and short, 
sharper flares appear with similar relative amplitude in the 
two energy bands. In order to measure any time delay, we 
calculated a CCF for each of the three RXTE orbits, without 
applying any de-trending procedure. The results are shown 
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Figure 2. Cross-correlations of the X-ray and IR light curves 
of GX 339-4 (positive lags mean IR lags the X-rays). A strong, 
nearly symmetric correlation is evident in all the three time in- 
tervals, corresponding to different RXTE orbits. In the inset we 
show a zoom of the peaks, showing the IR delay of ~100 ms with 
respect to the X-rays. The inset also shows a slight asymmetry 
toward positive delays. 



Figure 3. X-ray (2-15 keV) power spectrum of the second RXTE 
orbit (upper curve), together with the power spectrum of the 
simultaneous IR light curve (lower curve). The Poissonian noise 
has been subtracted from both spectra. The peak at ^6 Hz in the 
IR spectrum is instrumental. The high-frequency portion of the 
IR spectrum has yet un-modeled systematics, which however do 
not affect the results presented here. 



in Fig. [2] The strong correlation is confirmed. The CCF ap- 
pears highly symmetric and relatively stable over the three 
time intervals, with the change in amplitude simply reflect- 
ing the different variability amplitude in the light curves 
themselves. In the inset, we show a zoom on the peak of the 
CCF, which shows how the IR emission lags the X-rays by 
0.1 seconds, to which we associate an uncertainty of 30% 
(which includes systematics). 



4 DISCUSSION 

The main result of our work is the discovery of a strong 
correlation between the IR and the X-ray variability in GX 
339-4. The fact that the CCF is nearly symmetric and peaks 
at 100 ms rules out a reprocessing origin for the IR variabil- 
ity. If the IR radiation arose from reprocessing of X-rays by 
the outer disk, the short time delay would imply a highly in- 
clined disk. This would produce an highly asymmetric CCF, 
with a tail at long lags ( O'Brien et al.|[2002 1 . 

Additionally, power spectral analysis shows significant 
IR variability (at least 5% fractional rms, see Fig. [3| on 
timescales of ~200 ms or shorter, which sets an upper limit 
of ~ 6 x 10 9 cm to the radius of the IR-emitting region. From 
(5% of) the observed IR average flux of F ~ 1.5 x 10 -11 



erg s~ cm - , we derive a minimum brightness temperature 
of ~ 2.5 x 10 6 K. Optically thick thermal emission of the 
derived size and temperature would result in a 2-10 keV flux 
in excess of 10 -5 erg s _1 cm -2 , which is not observed in the 
data. These values represent very conservative estimates: a 
smaller region emitting the IR radiation would result in a 



higher brightness temperature, which would in turn result in 
a higher expected X-ray luminosity. With similar arguments 
we exclude thermal Bremsstrahlung emission. The existence 
of an IR lag is also inconsistent with the magnetic corona 
model ( jMerloni et al.|2000[ ), in which the same population of 
electrons produces the IR synchrotron emission and the X- 
ray Compton emission. We conclude that the most plausible 
origin for the observed IR variability is synchrotron emission 
from the inner jet j^] 

This result is a new, independent strong indication that 
jet synchrotron emission contributes significantly to the IR 
radiation in this source. This is the first time that hard- 
state, compact jet emission has been securely identified to 
vary on sub-second timescales in an XB, although variability 
on similar dynamical timescales toyn (i.e., scaled to mass) 
had been already observed in Active Galactic Nuclei (e.g. 
Scho deTet al.|2007| . These data thus represent a further step 
forward towards a full unification of the accretion/ejection 
process over a broad range of black-hole masses. 



4.1 Emitting regions and jet speed 

Our data strongly suggest that the variable IR emission 
comes from the jet, although we cannot conclude whether it 



2 This is confirmed by nearly-simultaneous optical and IR obser- 
vations, obtained while the source was in the same low-luminosity 
state. Those data (Lewis et al., in prep.) show a flat or inverted 
spectrum (inconsistent with thermal emission from a disc or com- 
panion star), and a long-timescale (>~minutes) variability stronger 
in IR than in optical. 
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is optically-thick or -thin synchrotron. The X-ray emission 
is usually interpreted as Comptonized radiation from ener- 
getic plasma in the very inner regions of the accretion flow, 
although the actual emitting region is still an open issue (ei- 
ther a corona or the base of the jet itself, for a discussion see 
|Markoff, Nowak fc Wilms|2005| |Maccarone|2005[ Malzac et 
al. 2009). Here we discuss the four possible scenarios. 

1 ) IR: optically thick - X-ray: inflow emission: the ob- 
served time delay between the IR and the X-ray variability 
gives an upper limit (given the unknown time for the ejection 
to take place) to the travel time of the variability - thus pre- 
sumably the matter - along the jet. Given a measure of the 
jet elongation we could estimate the jet speed. Such a mea- 
sure is not available for GX 339-4; however, a jet elongation 
measurement has been reported from 8.4 GHz observations 
of another BH XB, Cyg X-l ( |Stirling et al.|2001[). Within the 
standard model for compact jets ( jBlandford fc Konigl p979 
- BK79), the distance from the black hole of an emitting 
region in the jet is a function of the observed wavelength, 
the viewing angle, and the jet physical properties. 

Assuming that the main physical properties of the jet 
do not change, we can rescale the jet elongation measured 
at radio wavelengths in the BH Cyg X-l down to the IR 
wavelengths, obtaining a measurement of the distance of the 
IR-emitting region in the jet from the black hole in GX 339- 
Hcinz & Merloni (2004) showed that the small observed 



serve the flat spectrum - is largely untested (see e.g. Kaiser 
~ Pe'er fc Case lla 2009). However, we have chosen very 



2006 



scatter in the radio/X-ray fluxes relation (Gallo, Fender & 
Pooley 2003| ) implies that all stellar-mass BHs have very sim- 



ilar jet velocities, with a 3-sigma spread of A(/3F) < 1.6. (We 
note however that there is an increasing population of radio- 
quiet sources, whose nature is not yet understood (see e.g. 
|Gallo|[2007l |Casella fc P?erl[2009l Soleri et al. submitted.; 
and references therein). Their inclusion in the radio/X-ray 
flux correlation increases the scatter of the correlation itself, 
so the conclusions drawn by Heinz & Merloni ( |2004[ ) do not 
necessarily hold any longer.) 

We can use Eq. 28 of BK79 to obtain the GX 339-4 
jet elongation in IR, and thus the jet speed. We calculate 
the jet speed for a total of 10 sets of parameters, with each 
parameter randomly chosen within its measured permitted 



range. For Cyg X 1, we use a distance of 2.0±0.2 kpc (Gier- 
linski et al. 19 99]) a nd a source inclination over the 20°-70 c 



range ( polan|1992] |Ziolkowski|2005| ) . We use a projected jet 
elongation of Lj cy9 = 2.6 x 10 xDt pc cm (where Dk pc is 
the distance of Cyg X-l in kpc), to which we attribute a 50% 
uncertainty. For GX 339-4, we used the minimum distance 
of 6 kpc ( [Hynes et al.|2003| 2004) a nd a range of inclination 
angles 15°-60° ( [Cowley et al.||2002| |Miller et al.|2008[ ). 

We use the X-ray luminosity as a tracer of jet kinetic 
power, under the additional assu mption that the jet power 



follows the Pj oc L x relation (Fender, Gallo & Jonker 



20031. We estimate the X-ray (2-10 keV) fluxes of the two 



sources through spectral fitting of the RXTE data simulta- 
neous to the IR (GX 339-4, Fx = 1.4 x 10" 10 erg cm" 2 s" 1 ) 
and radio (Cyg X-l, Fx = 9 x 10~ 9 erg cm -2 s _1 ) obser- 
vations, and we attribute a 50% uncertainty to both values, 
to account for the uncertainties in the Pj — Lx scaling. 

We obtain a 3.3-ct lower limit of T > 2. It is important to 
remember that the reliability of this estimate depends on the 
key assumptions made in BK79 being true - the most impor- 
tant of these, that particle acceleration along the jet is con- 
tinuous and counterbalances energy losses in order to pre- 



conservative ranges for all parameters. The lowest Lorentz 
factors (2-3) are obtained only for small angles of Cyg X-l, 
for the minimum Cyg X-l jet elongation within the allowed 
range and for the minimum GX 339-4 distance. We con- 
clude that, if the IR synchrotron emission is optically-thick, 
these data suggest that the jets from accreting stellar-mass 
BHs are at least mildly relativistic, also in their common 
low/hard state. If, as is widely suggested, the jet speed cor- 
responds to the escape speed at the launch point, this might 
imply that the jet is launched from a region very close to 
the black hole itself. 

The obtained values of F might allow us further consid- 
erations. Assuming a random distribution of jet inclination 
angles 9 for the known population of stellar-mass BHs (i.e., 
Pe oc sin#), all with a Lorentz factor of 2, about 57% of 
them would appear de-boosted. This percentage increases 
up to 70% for r = 3, and up to ~77% for r = 4. Thus our 
jet speed estimate suggests that more than half of the jet 
luminosities measured from BHs in our Galaxy might have 
been underestimated. If this is true, there might be a pop- 
ulation of BH binaries with very high radio to X-ray ratios, 
which will be revealed in our and other galaxies by future 
radio all-sky monitors and large radio telescopes. The ob- 
tained estimate for F has several caveats, or at least large 
uncertainties. Future monitoring observations with the same 
technique will allow to refine this measure, studying the rel- 
ative dependency of this quantity with the varying accretion 
rate or total luminosity 

2) IR: optically thin - X-ray: inflow emission: this sce- 
nario is actually twofold: the IR optically thin emission could 
originate at the first shock in the jet, or further away af- 
ter cooling. In the first case, the observed time delay would 
measure the ejection and first acceleration timescale. In the 
second case, since we expect the spectral break to be around 
the IR wavelengths, we can approximate the IR radiation as 
if it was optically thick, and use Eq. 28 of BK79 as above. 
Thus, the lines of reasoning described in the previous sce- 
narios hold, as well as the lower limit to the jet speed. 

3) IR: optically thick - X-ray: jet emission: the reason- 
ing described in the first scenario, as well as the resulting 
lower limit for F, still hold, provided that the X-ray emission 
is not boosted (for a discussion about the evidence for this 
be the case, see |Heinz fc Merloni|2004| . 

4) IR: optically thin - X-ray: jet emission: the electron 
populations emitting at IR and X-ray wavelengths cannot 
be co-located, since we observe a delay, which here must 
represent a cooling time. If we assume that the X-ray emis- 
sion gives us a measure of the characteristic energy of the 
emitting electrons, we can place an upper limit on the mag- 
netic field intensity in the jet (see e.g. [Takahashi et al.|1996| 
for a use of this method in blazars). We obtain a unique 
solution for the system, with a magnetic field intensity of 
B~ 10 4 Gauss, an initial Lorentz factor of the electrons of 
7 ~ 7 x 10 3 , which becomes 7 ~ 50 after they cool down 



as to emit in IR. Electron re-acceleration (e.g. Jamil et al. 
| 2009[ ) would act as to counter balance the radiative losses, 
resulting in higher values for the magnetic field intensity. 
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4.2 Comparison with optical/X-ray CCFs 

The main difference between this IR/X-ray CCF and the 
published optical/X-ray CCFs (for this as well for other 
BHs) is the longer IR wavelength itself, which allows us to 
put strong constraints on the emission processes. In the hard 
state of BHs, the jet is expected to dominate the IR emis- 
sion ( Russell et al.|2006 see also 2 ). In fact, as we discussed 
in the previous section, these data allow us to rule out any 
thermal origin for the observed variability. To put similar 
constraints using optical data, a much higher luminosity (or 
faster variability) is needed (see e.g. Motch et al.|1982 ). 



Our CCF is very different from the one observed be- 



tween the optical and the X-ray variability (Gandhi et al. 
2008[ ). Namely, the sharp positive peak in our dataset is 
not preceded nor followed by any dip. This might be due to 
the additional presence of a large long-timescale variability, 
which buries the short-timescale structure of the IR CCF. In 
fact we did not apply any detrending procedure to our data, 
while this has been partially done for obtaining the optical 
CCF (Gandhi et al., in prep.). That this can be the reason 
for the observed differences seems to be suggested also by the 
fact that the IR CCF with the lowest amplitude among the 
three orbits (black curve in Fig. |5| appears steeper towards 
the positive time delays, similarly to the optical/X-ray CCF 
of this source. On the other hand, it might be that the the 
dips observed in the optical CCFs arise from a spectral com- 
ponent which does not contribute much at IR wavelengths. 
A detailed study of this will be presented in future work. 

The two CCFs differ also in the time delay measured at 
their peak: the optical/X-ray CCF peaks at ~ 150 ms, a 50- 
ms longer delay than IR (although we note that, depending 
on the uncertainty on the measured optical delay, the two 
delays might be marginally consistent within the errors). 
Given that the optical emission should come from closer to 
the black hole than the IR emission, this result is contrary to 
simple expectations. The two observations occured at simi- 
lar X-ray luminosities, casting doubt on jet power variation 
as a reason for the longer optical than IR delay. However, 
the optical data were acquired after the decay of a bright 
outburst, while our IR data were acquired at the end of the 
rise of a weak outburst, suggesting a hysteresis effect might 
play a role (jVadawale et al.| |2003 Fender, Bcllo ni fc Gallo| 
2004 Russell et al.|200"7 1. In particular, during/after the de- 
cay of an outburst the emission at a given wavelength might 
happen at larger distances from the black hole, because of 
the smaller amount of previously ejected matter that the 
jet encounters. Future simultaneous IR, optical and X-ray 
observations, at high-time resolution, will allow us to unveil 
some of these important unknowns, thus allowing us to test 
some of the physical assumptions which are now behind this 
method. 
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